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Abstract. The characterization of the growing number of newly discovered exo- 
planets — nature, internal structure, formation and evolution — strongly relies on the 
properties of their host-star, i.e. its mass, radius and age. These latter can be inferred 
from stellar evolution models constrained by the observed global parameters of the 
host-star — effective temperature, photospheric chemical composition, surface gravity 
and/or luminosity — and by its mean density inferred from the transit analysis. Addi- 
tional constraints for the models can be provided by asteroseismic observations of the 
host-star The precision and accuracy on the age, mass and radius not only depend on 
the quality and number of available observations of the host-star but also on our ability 
to model it properly. Stellar models are still based on a number of approximations, they 
rely on physical inputs and data that can be uncertain and do not treat correctly all the 
physical processes that can be at work inside a star. We focus here on the determi- 
nation of the age of HD 17156, an oscillating star hosting an exoplanet. We examine 
the dispersion of the age values obtained by different methods — empirical or model- 
dependent — and the different sources of errors — observational or theoretical — that 
intervene in the age determination based on stellar models. 



1. Introduction 

With the advent of ultra high precision photometry as performed on board by the spatial 
missions CoRoT and Kepler and the HST, it is now becoming possible to detect with 
the same device the transit of an exoplanet in front of its host-star and the host-star 
oscillations. The analysis of the light curve, which carries the signature of the transiting 
object, provides the normalized separation a/R^, where a is the separation between the 
two bodies and the radius of the star. This quantity is a measure of the mean host- 
star density (pi,), through the third Kepler law. The light curve analysis also provides 
the planetary radius Rp in units of 7?^, the impact parameter a cos///?* and the period 
of the motion. Ground-based high resolution spectroscopy can provide the star radial 
velocity curve and in turn the planetary mass Mp in units of Mi,. It is then necessary to 
estimate the mass and radius of the host-star to get the absolute mass and radius of the 
exoplanet and classify it as a gaseous giant body or as a telluric solid one. Furthermore 
to get insights on the planet internal structure, formation and ev olution, we need to 
know not only the radius and mass of the host-star but also its age ( Havel et"al]|201 1 ). 



Stellar ages cannot be determined by direct measurem ents but can only be esti 



mated or inferred. As nicely reviewed by ISoderblomI (1201(1) they are many methods to 
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estimate the age of a star according to its mass, age and configuration — single star or 
star belonging to a group. Here we focus on the determination of the age of HR 17156, 
a pulsating star th a t host s an exoplanet. This star w as observed and modeled recently 
by iNutzman et alj ( 2011 ) and iGilliland et alj (l201lb who estimated its age from com- 
parison with stellar models. In Section [2] we present the observational data available 
for HR 17156. In Section[3]we present the empirical methods available to age-date the 
star and we revise and compare the resulting ages. In Section 0] we discuss the ages 
obtained through isochrone placement in the H-R diagram. Then we proceed with a 
detailed modeling of the star. We concentrate on quantifying the sources of inaccuracy 
affecting the age inferred from stellar models which result from the choice of the in- 
put physics and parameters. In Section [5] we synthesize the results and emphasize the 
progress allowed by asteroseismology. 



2. HD 17156: Observational Data and Constraints for the Modeling 

HD 17156 (HIP 13192) is a si ngle G5V s tar, located at w75 pc according to its Hippar- 
cos parallax (13.33±0.72 mas. Ivan Leeuwe n 2007). We consider the observational data 
of HD 17156 listed in Table [T] We adopt the sam e effective temperatur e Tgff, surface 
gravity log^, photospheric metallicity [Fe/H] than'Gilliland et all ( 20 lib . We have de- 



rived the luminosity L/Lq from the para llax and the V-magnitude y =8.17±0.01, with 
the bolometric correction BC=-0.016 of lVandenBerg & ClemI (l2003h . 

Ultra high precision photometry performed on HST has provided a light-curve car- 
rying th e signature of an exopl anet and of stellar oscillations. The transit has been anal- 
ysed by iNutzman et all (l201lh who derived a mea s ure of the mean density of the star 
which we express here as M^^^/?~ MGilliland et alj ( 20TTh analyzed the observed oscil- 



lation spectrum, a typical p-mode solar-like one and i dentified 8 lo w-degree p-modes 
(see their Table 1). In the asymptotic approximation ( Tassoul|[l980 ) the frequency of 



a mode of radial order n and angular degree writes Vnj=Av{n+^£+e^ -£(£+l)Do. 
The quantity e is sensitive to surface physics but weakly sensitive to the order and de- 
gree of the mode. The difference in frequency between modes of consecutive orders 
and same degree is approximately constant and given by Av— v„4-i^— y^/=Avf while the 
difference in frequency between modes of consecutive orders and degrees differing by 
two units is 5y^=v„/-v„_i/+2=4(/' -i- 6)Do. The differences Ay and dv are called the 
large and small separations, respectively. The large frequency separation Ay is a mea- 
sure of the inverse of the sound travel time across a stellar diameter and is related to 
the mean stellar density. The quantity Dq is sensitive to the sound speed gradient in the 
inner regions which changes with evolution and th erefore with age . The values of Ayo, 
Do and e derived from the fit of the frequencies by IGilliland et al.l are listed in Table [T] 



Table 1. Observational constraints for HD 17156 



Teff 




[Fe/H] 


L 




Avo 


Do 


e 


[K] 


[dex] 


[dex] 


[Lo] 




[juHz] 


[//Hz] 


[//Hz] 


6082 


4.31 


0.24 


2.45 


0.718 


83.44 


0.90 


1.15 


60 


0.04 


0.03 


0.28 


0.007 


0.15 


0.19 


0.04 
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3. Age of HD 17156 from Empirical Methods 



We discuss briefly here the empirical methods that can be used to estimate the age of 
HD 17156 and we report the results in Figure [T] (right figure, second column). 

Activity. The chromospheric activity of solar-type dwarfs is anti-correlated with their 
age. Empirical relations allow to rely the Call H & K emission index R'^^ = LnK/^boi 
to age (see e.g.'Mamajek & Hillenbrandll2008l . for a recent calibration). For HD 17156, 
iRscher et al. (2007) measured log 7?^ ^ = -5.04 which smallness indicates very low 
chromospheric activity. Using |Mamaiek & Hillenbran J'j?^ -age relation, we find an 
age of 7.4±4.3 Gyr. Also from the Mamajek & Hillenbrand relation between the frac- 
tional X-ray emission R'^ = Lx/L\,oi and age, with Lx<28.7 (.Kashyap et al..i2008i) . we 
estimate a lower age limit of 1.5 Gyr for the star. 



Gyrochronology. When they evolve, solar-type stars lose angular momentum via 
magnetic breaking due to thei r mass l oss. It leads to a decrease of their rotation rate, 
first quantified by ISkumanichI d 19721) . 



BamesI (120071) has proposed a new method 



gyrochronology — to derive the age of solar-type s tars from the empirical relation 
linking their rotation pe r iod, c olor and age. From iBamesl 's relation as revised by 
Mamaiek & Hillenbrandl ( 20081) and the rotation period Prot-12.8 days measured by 
Fischer et al. (2007), we derive an age of 1.4±0.3 Gyr for HD 17156. 



Photospheric Lithium Abundance. At the surface of low mass stars, the lithium 
abundance can be depleted when the convective zone reaches the not so deep regions 
where Li is destroyed by nuclear reactions at 7*2.5 x 10^ K or if mixing processes 
carry Li from the basis of the convective zone to the nuclear burning region. A relation 
between the Li -abundance, effe ctive temperature and age is observed (but not fully un- 
derstood). Ba rbieri et al.l ( 20091) measured the Li abundance of HD 17156 and inferred 
a lower age limit of 2 Gyr on the bas is of Li abundance curves as a function of Teg 
published by ISestito & RandichI (|2005|) for clusters of different ages. 



4. Age of HD 17156 from Model-dependent Methods 

The age and mass of stars can be inferred from stellar evolution models. The mass 
and age (i.e. the inputs of models) have to be adjusted in order to get a model that 
fits the observed parameters of the host-star. The more observational constraints, the 
better determined mass and age. We consider here two methods to age-date HD 17156: 
datation based on the placement of the star on pre-calculated isochrone grids — which 
is a thoroughly used method — and datation by means of a complete modeling of the 
star. 

4.1. Isochrone placement 

This method consists in placing the star in a H-R diagram and interpreting its posi- 
tion by means of a grid of pre-calculated theoretical isochrones or evolutionary tracks. 
In addition to the fact that the star position is defined modulo the observational error 
bars on luminosity, effective temperature and metallicity, this inversion method suf- 
fers from several uncertainties. First, isochrones grids are derived from stellar models 
that are based on a number of assumptions/simplifications and that include uncertain or 
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badly known input physics. Second, observations cannot be compared directly to the- 
oretical isochrones. In particular model atmospheres, also based on assumptions and 
approximated physical descriptions, have to be used to derive the effective temperature, 
bolometric magnitude and chemical abundances of stars from spectra or colors. Fur- 
thermore different inversion techniques can be used to extract the stellar age and mass 
from theoretical isochrones but they have to deal with problems of degeneracy in some 
regions of the H-R diagram, f or instance when th e star is close to the zero age main 
sequence or at turn-off (see e.g.' Pont & Eyeill2004t) . 

As for HP 17 156, we derive an age of f*=3.2±0.2 Gyr from Padova isochrones 
(iGirardi et alJl2002h using a Baye s ian inversion techniqu e similar to the one designed 
by |j0rg ensen & LindegrenI ( 2005 ). iNutzman et al.l ( 20111) estimated an age of S.SS^q ^^' 
Gyr by isochrone placement but incorporating the constraint on the mean-density from 
the transit. To summarize we show in Fig. \T\ (first column), several ages values of 
HD 17156 derived from isochrones inversion and picked up in the literature. Notewor- 
thy, the ages obtained cover a large range, 2.3-7.2 Gyr, due to the different isochrones 
grids and inversion methods used. 
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Figure 1. Age (left figure) and mass (right) estimates for HD 17156. Column 
"HR" shows estimates based on (i) inversion of isochrones by us (full grey circle) or 
reported from the literature (open circles) and (ii) complete modeling of the star 
with the H-R diagram constraints only (model AO, full black). Column "Emp" 
displays empirical ages. Column "Transit" shows values from complete modeling 
based on the H-R diagram and transit constraints (model BO, full black circle). Col- 
umn "Seismo" are values from complete modeling based on global and seismic con- 
straints and different sets of input physic s (model CO, ful l black circl e, models CI -7 
open c ircles). Squares are values from iNutzman et al.l (1201 lb and iGilliland et al.l 
(l20TTh 



4.2. Complete Modeling 

We h ave modeled the star with the stellar evolution code Cesam2k dMorel & Lebreton 



20081) and calculated the oscillations frequencies with the LOSC pulsation code (IScuflaire et al 
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2008b . We have used the Levenberg-M arquardt minimization method in the way de- 
scribed by 'Mi glio & MontalbanI (1200 5*) to adjust the unknown parameters of the mod- 
eling so that the model of HD 17156 fits at best the observations, within the error bars. 
Our models are based on the input physics and parameters described below. If enough 
observational constraints are available, the free parameters can be adjusted in the mod- 
eling process, otherwise they have to be fixed — ^more or less arbitrarily by the modeler 
(see below). 

• Microsc opic Physics: The referenc e models are based on the OPAL05 equation 



of sta te (iRogers & Nayfonovll2002l) and on OPAL96 opacities (Iglesias & Rogers 



1996) complemented at low temperatures by W ICHITA tables (.Ferguson et al 



20051) . We use the NACRE nuclear reaction rates (lAngulo et al.ll 19991) except for 



the ^^N{p, 7)^^0 reaction where we adopt the revised LUNA rate (For micola et al 



12004 ). Models take into account microscopic diffusion of helium and heavy ele- 
ments — including gravitational settling, thermal and concentration diffu s ion bu t 
no radiative levitation — following the formalism of iMichaud & Proffittl(ll993h . 
In alternate models we consider the OPALOl equation of state, the NACRE nuclear 
reaction rate for the ^^N{p, y)^^0 reaction and assume no diffusion. 



Macroscopic Physics: We use the COM convection theory of iCanuto et al.l (1 1 9961) 
with a free mixing-length parameter aconv We consider the MLT theory dBohm-Vitense 
[1958) as an alternative. We include convective core overshooting and assume that 
the overshooting zone is adiabatic and fully mixed. As a reference, we set the 
overshooting distance to be ^ov-Q^ov x niin(/?cc> ^p) where Oov, Rcc and Hp are 
the overshooting parameter, the radius of the convective core and the pressure 
scale height respectively. In alternate models, we adopt the Roxburgh (199_^ 
prescription, in which overshooting extends on a fraction of the mass of the con- 
vective core Mcc, the mass of the mixed core being expressed as Mov^aov x ^cc- 

Atmospheric Boundary Condition: The reference models are based on grey model 
atmospheres with the classi cal Eddington T-r law but we also investigate models 
based on the Kurucz (1 1 9931) T-T law. 



• Solar mixture: We adopt the revised AGSS09 solar m ixture ( Asplund et al.l 20091) 
as the reference but we consider the GN93 mixture dGrevesse & Noelslll993h in 
an alternate model. 

As described below we have considered different situations where k unknowns of 
the model are adjusted to fit k observational constraints. 



Age from Global Parameters. In the modeling case AO we suppose that solely the 
global parameters are constrained by observation (rgff, present [Fe/H], L) and we seek 
which mass, age and initial metal to hydrogen ratio Zq/Xq are required for the model 
to satisfy these constraints. Since this gives 3 unknowns for 3 observed pai^ameters we 
have to make assumptions on the other inputs of the models, mainly the initial helium 
abundance Yq, the mixing-length and overshooting parameters. We assume that Yq 
can be derived from the helium to metal enrichment r atio (Yq - Yp) /Z=AY / AZ where 
yp=0.245 is the primordial He abundance (e.g. Pe imbert et al.ll2007[) and AF/AZ«2 is 
taken from a solar model calibration in luminosity and radius. We take (2conv=0.768 
from the solar model calibration and use a moderate overshooting parameter aov=0.15. 
The input physics are the reference ones (hereafter denoted as REF). 
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Age from Global Parameters and Transit. In the modeling case BO we add the 
constraint on MiR^^ — derived from the planet transit — to the global parameters con- 
sidered in case AO and we adjust 4 unknowns of the models (age, mass, metallicity and 
CKov) with aconv and Yq fixed as in Case AO. The input physics are the reference ones. 

Age from Global Parameters, Transit and Asteroseismology. In the different cases 
CO-7, we add the seismic constraints Avq, Dq to the modeling and we adjust 6 model 
unknowns (age, mass, Yq and Zq/Xq, Ocom, Q^ov)- We consider several alternatives for 
input physics: Case CO is for reference inputs. Case CI uses the the NACRE reaction 
rate for the ^'^N{p, y)'^0 rate. Case C2 uses the MLT convection formalism. Case C3 is 
without diffusion. Case C4 uses the OPALOl equation of state. Case C5 uses Kurucz' 
model atmospheres as boundary condition, the MLT formalism and the GN93 mixture. 
Case C6 models overshooting as being a fraction of the convective core mass. Case C7 
is based on the GN93 mixture. 



Table 2. Parameters of HD 17156 inferred from complete modeling 



Case Inputs age M R Yq Zq/Xq ffconv ^ov 

[Gyr] [Mo] [R^] 



AO 


REF 


4.33 


1 


22 


1 


41 


0.299 


0.0376 


0.77 


0.15 


10-4 


BO 


REF 


4.01 


1 


26 


1 


50 


0.298 


0.0376 


0.77 


0.08 


1.3 


CO 


REF 


3.91 


1 


26 


1 


50 


0.297 


0.0383 


0.72 


0.10 


1.7 


CI 


NACRE 


4.06 


1 


25 


1 


49 


0.298 


0.0370 


0.73 


0.16 


1.6 


C2 


MLT 


3.67 


1 


25 


1 


50 


0.298 


0.0374 


1.76 


0.15 


2.1 


C3 


no diffusion 


4.60 


1 


25 


1 


49 


0.283 


0.0308 


0.69 


0.26 


3.7 


C4 


OPALOl 


3.87 


1 


26 


1 


50 


0.295 


0.0382 


0.78 


0.10 


3.2 


C5 


Kurucz 


4.18 


1 


23 


1 


49 


0.300 


0.0466 


1.95 


0.13 


6.6 


C6 




3.92 


1 


26 


1 


49 


0.297 


0.0380 


0.73 


1.18 


1.5 


C7 


GN93 


3.63 


1 


28 


1 


50 


0.300 


0.0498 


0.73 


0.05 


1.6 



Table|2]lists the values of the parameters of HD 17156 estimated through complete 
modeling and the corresponding minimum value. The ages and masses obtained in 
the different cases are displayed in Fig.H] In left and central Figures |2] we show that the 
reference model CO satisfies quite well the observed position of HD 17156 in the H-R 

diagram and in the M^R'^-T^s plane. The echelle- diagram is displayed i n Figure |2] 
(right). It shows that the 8 frequencies identified by 'Gillilan d et all ( 201lh are well- 
fitted by the CO model for the corresponding value of their degree. 



5. Synthesis and Conclusions 

We have investigated several methods to age-date the star HD 17156. We have found 
that the age estimated through empirical techniques is very ill-defined for this star which 
shows quite low activity. On the other hand, the age f*=3.2±0.2 Gyr derived from the 
placement of the star on a set of pre-calculated isochrones — given its luminosity, Teff 
and [Fe/H] — has an internal error of about 10 per cent. However quite a large scatter 
(w3.7 Gyr) is found when comparing the ages found in the literature because they are 
based on different observational parameters, inversion methods and isochrones grids. 
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Figure 2. Fit of the observational constraints by the Case CO reference model. 
Left: Evolutionary track in the H-R diagram with respectively the Icr and 2cr errors 
bars on luminosity and temperature. The dot locates the CO-model, Centre: Similar 
to the previous one but for the M"/? '-Feff plane. Right: Echelle diagram with on 
ordinates the frequencies of the modes and on abscissas the same frequencies modulo 
the large separation. Open symbols correspond to the frequen cies predicted by mode l 
CO while full symbols correspond to the 8 modes identified by lGilliland et al.l (l201 Ih . 
Circles, triangles and squares are for degrees {-Q, 1 , 2 respectively. 




'3.81 3.8 3.79 3.78 3.77 3.76 3.8 3.79 3.78 3.77 20 40 60 80 

log T„ log T,„ " modulo Ak [rfia] 



The complete modeling of the star presented here provides different ages depend- 
ing on the number of observational constraints considered in the modeling. In Case AO, 
when only 3 constraints are considered (luminosity, T^s, Fe/H), assumptions have to 
be made on the values of important model parameters (initial helium abundance, over- 
shooting and mixing length parameters). Interestingly these latter are determined by 
the calibration in Cases CO-7 when seismic and transit constraints are considered. The 
difference in age between Case CO (3.91 Gyr) and Case AO is of 0.4 Gyr, i.e. about 
10 per cent. Furthermore we have found a scatter in age of about 25 per cents (1 Gyr) 
when using different input physics in the models (Cases Cl-7). This should be kept in 
mind when giving the age of stars from a given set of stellar models 

By-products of the calibration are the mass M^^ and radius of the star. As 
shown in Table |21 the differences in and resulting from variations of the inputs 
of the models are of less than 3 and 1 per cents respectively. The calibration also 
yields a rather low valu e of t he overshooting parameter Qfov=0.10 which according to 
Roxburgh's formalism ( 19921) corresponds to an extension of mixing on a fraction of 



the mass of the convective core M^c of 18 per cents. This will be investigated with more 
details in a forthcoming paper. 

Acknowledgments. I warmly thank Josefina Montalban, Arlette Noels and Math- 
ieu Havel for valuable advice, documents and fruitful discussions. I am very grateful to 
John Leibacher who presented this talk for me in Hakone. 
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